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Abstract. We present supplementary information on the recent indirect measurement of the Casimir pres- 
sure between two parallel plates using a micromachined oscillator. The equivalent pressure between the 
plates is obtained by means of the proximity force approximation after measuring the force gradient be- 
tween a gold coated sphere and a gold coated plate. The data are compared with a new theoretical approach 
to the thermal Casimir force based on the use of the Lifshitz formula, combined with a generalized plasma- 
like dielectric permittivity which takes into account interband transitions of core electrons. The theoretical 
Casimir pressures calculated using the new approach are compared with those computed in the framework 
of the previously used impedance approach and also with the Drude model approach. The latter is shown to 
be excluded by the data at a 99.9% confidence level within a wide separation range from 210 to 620 nm. The 
level of agreement between the data and theoretical approaches based on the generalized plasma model, or 
the Leontovich surface impedance, is used to set stronger constraints on the Yukawa forces predicted from 
the exchange of light elementary particles and/or extra-dimensional physics. The resulting constraints are 
the strongest in the interaction region from 20 to 86 nm with a largest improvement by a factor of 4.4 at 
26 nm. 



1 Introduction 

It is well known that there is little or no experimental 
confirmation for many predictions of unified field theo- 
ries, supersymmetry, supergravity, or string theory. Direct 
experimental tests for many of these predictions require 
accelerators of very high energies which will be not avail- 
able in the foreseeable future. For this reason any non- 
accelerator tests of the predictions of new physics beyond 
the standard model attract the serious attention of both 
experimentalists and theorists. 

One of the most intriguing predictions made by many 
extensions of the standard model is the existence of light 
and massless elementary particles which arise as a result 
of some spontaneously (or weakly dynamically) broken 
symmetry. Beams of such particles can penetrate through 
thick matter with practically no interaction. This makes 
it difficult to investigate these particles using the usual 
laboratory setups of elementary particle physics. There is, 
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however, an alternative way to investigate light elemen- 
tary particles and their interactions by using table-top 
laboratory experiments. These experiments utilize the fact 
that the exchange of such particles between atoms belong- 
ing to two different macrobodies can generate a new long- 
range force in addition to the commonly known electro- 
magnetic and gravitational interactions. For example, the 
exchange of predicted light bosons, such as scalar axions, 
graviphotons, hyperphotons, dilatons and moduli among 
others (see, e.g. 1,2,3,4,5J) generates a Yukawa potential. 
The simultaneous exchange of two photons, two massless 
scalars or massless pseudoscalars, and the exchange of a 
massless axion or a massless neutrino-antineutrino pair 
leads to power-law interactions with different powers [6] 
[TUSM HUlflT] . Coincidently, a Yukawa correction to New- 
tonian gravity is predicted in extra-dimensional physics 
with compact extra dimensions and low energy compacti- 
fication scale of order of 1 TeV [T2Tfl~3T[L4l[l~5] . Furthermore, 
some brane theories contain exactly the standard model 
at low energy [16j . For models of non-compact but warped 
extra dimensions, power-law corrections to the Newtonian 
gravitational law have been predicted [17] . The cosmolog- 
ical constant generated in such models may be of the cor- 
rect order of magnitude as suggested by observations [18] . 
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Direct experimental signatures of strings and branes are 
discussed in [T9] . 

Experimental constraints on hypothetical long-range 
interactions arising from both light elementary particles 
and large extra dimensions can be obtained from precise 
force measurements between macrobodies. For electrically 
neutral test bodies the dominant background force at sep- 
arations greater than 10 -5 m is the gravity. At shorter sep- 
arations the dominant forces are the van der Waals and 
Casimir forces caused by fluctuations of the electromag- 
netic field [20 . During the past few years a number of new 
experiments have been performed to measure small forces 
between macrobodies and to obtain stronger constraints 
on hypothetical long-range interactions (which are also re- 
ferred to as the "fifth force" pQ). Thus, in sub-millimeter 
gravity experiments stronger constraints on Yukawa cor- 
rections to the Newtonian gravitational force for ranges 
~ lfr 4 m and ~ 10" 5 m have been obtained [2l1l2"2"Il2"3TI2H 
25, 26J . In a series of experiments measuring the Casimir 
force between gold coated test bodies the constraints on 
Yukawa-type interactions in a sub-micrometer range have 
been strengthened up to 10 4 times [27 | [2^ p9l[30l l3T 1 [32 l l33 l 

[SUMMIT]. 

This paper exploits the results of the most precise re- 
cent determination of the Casimir pressure between two 
parallel gold coated plates using a micromechanical tor- 
sional oscillator. This is the third in a series of experiments 
using a micromechanical oscillator for precise Casimir force 
measurements. Results of the first two experiments were 
published in [36,37 (previously a similar technique was 
used to demonstrate the actuation of a micromechanical 
device by the Casimir force [38]). A brief discussion of the 
results of the third experiment, and a description of the 
main improvements, as compared with the previous two 
experiments, is contained in [39j . Here we present addi- 
tional experimental details related to the experiment [39) 
which were not discussed in the first publication, includ- 
ing the resistivity measurements and tests of the linear- 
ity of the oscillator used. The focus of this paper is a 
comparison of the experimental data with a recently pro- 
posed new theoretical approach to the thermal Casimir 
force [40] which is applicable to all experiments regardless 
of the separation between the interacting bodies. Within 
this framework, we first present a precise fit of the tabu- 
lated optical data [41] for the imaginary part of the dielec- 
tric permittivity of gold within a wide frequency region. 
The fit is obtained using a set of six oscillators represent- 
ing interband transitions in gold. We then compare this 
fit with a previously known fit based on DESY data [42J, 
[43] . Our theoretical approach based on the Lifshitz for- 
mula is found to be in very good agreement with the mea- 
sured results. The same measured results are also com- 
pared with an alternative approach to the theory of the 
thermal Casimir force [44], which approach is found to 
be excluded by our measurements at a confidence level of 
99.9%. The level of agreement between our theory and ex- 
perimental data is used to set constraints on Yukawa-type 
corrections to Newtonian gravity originating from the ex- 
change of light hypothetical elementary particles and/or 



extra-dimensional physics. The resulting constraints are 
several times stronger than those derived from previous 
experiments. We also reanalyze constraints following |35j 
from experiment [30] (in [35] the confidence level of our 
results was not determined). As a consequence, the inter- 
action region where the constraints from the present ex- 
periment are the strongest is widened. Special attention is 
paid to minor deviations between experiment and theory 
at the shortest separations. Although these deviations are 
inside the error bars and thus not statistically meaningful, 
we present an analysis of various explanations for them. 

The plan of this paper is as follows: in Sec. 2 we present 
a brief description of the experimental setup and mea- 
surement results with an emphasis on novel aspects not 
described previously in [36 |f3~7ll39] . Sec. 3 is devoted to the 
comparison of experimental data with different theoret- 
ical approaches including our new approach in [40]. The 
new precise oscillator fit of the optical data for gold is also 
presented here. Sec. 4 contains constraints on hypothetical 
Yukawa interactions following from the level of agreement 
of data with theory, and includes a comparison with con- 
straints obtained from earlier experiments. In Sec. 5 we 
present our conclusions and discussion. 

2 Experimental setup and measurement 
results 

One component of our setup is an Au-coated sapphire 
sphere attached to an optical fiber. The thickness of the 
Au coating on the sphere is Z\g = 180 nm, and the radius 
of the coated sphere is R = 151.3 ± 0.2 /xm. The sphere is 
placed at a separation z above a micromachined oscillator 
consisting of a heavily doped, Au-coated polysilicon plate 

(the thickness of the coating is Ag = 210 nm) suspended 
at two opposite points by serpentine springs. This plate 
can rotate under the influence of the Casimir force F{z) 
acting between the sphere and the plate. The rotation 
angle is measured by the change of the capacitance be- 
tween the plate and two independently contacted polysil- 
icon electrodes located under it (details of the setup are 
described in [36,37]). The micromachined oscillator and 
the sphere with a fiber were mounted inside a can with 
magnetic damping vibration isolation, where a pressure 
below 10 -4 torr was maintained. 

In this experiment a dynamic measurement mode was 
employed. For this purpose the vertical separation be- 
tween the sphere and the plate was varied harmonically, 
z(t) = z + A z cos(w r i), where u> r is the resonant angular 
frequency of the oscillator in the presence of the sphere. 
The magnitude of the amplitude A z « 2 nm was chosen in 
such a way that the oscillator exhibited a linear response. 
In the presence of the Casimir force F(z) the resonant 
frequency uj r is shifted relative to the natural angular fre- 
quency of the oscillator ujq = 2n x (713.25 ± 0.02) Hz de- 
termined in the absence of the sphere. In the linear regime 
this shift can be found using [3H l [37 ] [38 l [39] 

Ji _ 2 T 1 6 2 9F(z)l 
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where b is the lever arm between the axis of plate rotation 
and the projection on the plate of the closest point of the 
sphere, / is the moment of inertia of the oscillator, and 
b 2 /I = (1.2432 ± 0.0005) fig- 1 . 

The actual measured quantity in this experiment is the 
change of the resonant frequency of the oscillator, 
under the influence of the Casimir force F(z) acting be- 
tween the sphere and the plate. Using ([I}, the experimen- 
tal data for uj r — ojq obtained at different separation dis- 
tances can be transformed into dF(z)/dz. It is less useful, 
however, to recover the force F(z) between a sphere and 
a plate using the force gradient. A better avenue is given 
by using the proximity force approximation (PFA) 45,46, 

B3 



F(z) = 2irRE{z), 



(2) 



where E{z) is the Casimir energy per unit area between 
two infinitely large parallel plates composed of the same 
materials as the sphere and the plate. Differentiating with 
respect to z and taking into account that the Casimir pres- 
sure between the two parallel plates is 



P{z) 

one arrives to the expression 



dE{z) 



dz ' 



P{z) 



1 dF{z) 
2-kR dz 



(3) 



(4) 



From equations |T]) and (H|) one can immediately convert 
the experimental data into data for the Casimir pressure 
between two parallel plates. This is in fact the so-called 
indirect measurement 48J of the pressure. Note that in 
[4"9"] , where the configuration of two parallel plates was 
actually used in the experimental setup, the directly mea- 
sured quantity was also the frequency shift due to the 
Casimir pressure proportional to dP(z)/dz. The pressure 
P{z) was then recovered using the data for its derivative. 

The calibration of absolute separations between the 
plate and the sphere was performed by the application of 
voltages in a manner analogous to that reported in [36l 
[37]. The use of a two-color fiber interferometer [50] and 
a « 7% improvement in vibration noise yielded an error 
of only 0.2 nm in a distance z meas between the end of the 
fiber and the stationary reference. As a result, for every 
repetition of the Casimir pressure measurement we were 
able to reposition our sample to within Az meas = 0.2 nm. 
Finally the absolute separations z between the sphere and 
the plate were measured with an absolute error Az = 
0.6 nm determined at 95% confidence [37] . 

The indirect measurements of the Casimir pressure 
Pjiz-i) were repeated at practically the same separations 
%i (1 < i < 293) 33 times (1 < j < 33). The mean values 
of the experimental Casimir pressure 



33 



(5) 



are plotted in Fig. 1 as a function of separation over the 
entire measurement range from z\ — 162.03 nm to Z293 = 



745.98 nm. As an example, a few mean Casimir pressures 
P(zi) at different separations are presented in column (a) 
of Table 1. In this measurement the random experimental 
error is much smaller than the systematic error. Specifi- 
cally, using Student's i-distribution [5T] with a number of 
degrees of freedom / = 32, and choosing ft = 0.95 confi- 
dence, we obtain p = (1 + l3)/2 = 0.975, and t p (f) = 2.0. 
This leads to the random experimental error 



A^P^izi) = s( Zi )t p (f), 



(6) 



where s(z) is the variance of the mean for the pressure 



s 2 ( Zl ) = 



1 



33 

1056 ^ 



[P 3 {z t )-P{ Zl )\ 



(7) 



The random error in ([6]) reaches a maximum value equal 
to 0.46 mPa at z = 162 nm, decreases to 0.11 mPa at z = 
300 nm, and maintains this value up to z — 746 nm. 

The systematic error of the pressure measurements in 
this experiment is determined by the errors in the mea- 
surements of the resonance frequency, radius of the sphere 
(these errors were indicated above) , and also by the error 
in using the PFA. Until 2006 the latter was not known 
with certainty but estimated to be of order z/R on the 
basis of dimensional considerations [47]. Recently, how- 
ever, quantitative results on the accuracy of PFA were 
obtained theoretically for the configuration of a cylinder 
above a plate 52,53J (the electromagnetic Casimir effect), 
and for a sphere above a plate ,54,155] (the scalar Casimir 
effect). In addition the validity of the PFA was estab- 
lished experimentally [56] for a sphere above a plate. In 
all cases at small separations the error in using the PFA 
was shown to be less than z/R. However, in our conser- 
vative error analysis we estimate this error with a safety 
margin as z/R. By combining all the above systematic er- 
rors at 95% confidence using the statistical rules described 
in [37], we obtain a systematic error equal to 2.12mPa at 
z = 162 nm. The systematic error decreases to 0.44 mPa at 
z = 300 nm, and then to 0.31 mPa at z = 746 nm. Finally 



P(Pa) 




200 300 400 500 600 700 

z (nm) 

Fig. 1. Values of the mean Casimir pressure between two Au- 
coated plates as a function of separation. 
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Table 1. Magnitudes of the mean experimental Casimir pres- 
sures P (column a) at different separations z compared with 
the magnitudes of the theoretical pressures P th computed us- 
ing the generalized plasma model approach (column b), the 
Leontovich surface impedance approach (column c), the Drude 
model approach (column d), and with the half- width, S, of the 
95% confidence interval for P th — P (column e). All pressures 
are given in mPa. 



Z (nm) 


a 


b 


c 


d 


e 


162 


1108.4 


1098.4 


1094.2 


1076.2 


21.2 


166 


1012.7 


1007.1 


1002.7 


985.40 


19.0 


170 


926.85 


923.71 


919.56 


902.96 


17.1 


180 


751.19 


750.58 


747.06 


732.14 


13.3 


190 


616.00 


616.71 


613.70 


600.28 


10.5 


200 


510.50 


511.26 


508.70 


496.62 


8.40 


250 


225.16 


225.71 


224.45 


217.11 


3.30 


300 


114.82 


114.87 


114.18 


109.48 


1.63 


350 


64.634 


64.574 


64.176 


61.004 


0.98 


400 


39.198 


39.096 


38.850 


36.617 


0.69 


450 


25.155 


25.034 


24.874 


23.247 


0.54 


500 


16.822 


16.785 


16.678 


15.456 


0.47 


550 


11.678 


11.669 


11.595 


10.654 


0.42 


600 


8.410 


8.365 


8.312 


7.573 


0.39 


650 


6.216 


6.151 


6.113 


5.522 


0.38 


700 


4.730 


4.626 


4.598 


4.118 


0.36 


746 


3.614 


3.620 


5.598 


3.198 


0.35 



we combine the resulting random and systematic errors at 
a 95% confidence to arrive at the total experimental error, 
Z\ tot P cxp (z), approximately equal to the systematic error 
at all separations considered. Detailed information on the 
statistical methods used in our error analysis can be found 
in [3"Tl[4"B] . As a result, the total relative experimental er- 
ror AP ex P(z)/|F(z)| varies from 0.19% at z = 162 nm, to 
0.9% at z = 400 nm, and to 9.0% at z = 746 nm. Hence 
this is the most precise experiment on the Casimir effect 
performed up to date. 

Several additional measurements and tests were per- 
formed in order to compare the experimental data with 
theory in a conclusive manner. In order to include the ef- 
fects of surface roughness in theoretical computations of 
the Casimir pressure we have investigated the topography 
of the metallic coatings both on the plate (p) and on the 
sphere (s) using an AFM probe in tapping mode. All AFM 
scans were squares with sizes ranging from 0.5x0.5 /jm to 
10 x 10 /im. The information obtained was indistinguish- 
able. In the case of a sphere the surface curvature was 
taken into account. For this purpose the image was pla- 
narized, and then the roughness analysis performed. For 
a typical scan of 5 x 5 /im the effect of curvature is about 
40 nm. From AFM images of the surfaces, the fraction of 
each surface area vf with height h^'^ was determined. 
It was found that for the sphere (1 < i < = 106) 

h\ s ^ varies from to 10.94 nm and for the plate (1 < 
i < ijfW = 85) varies from to 18.35 nm. Here, the 
highest peaks on the sphere and on the plate are almost 



of the same height as in the previous experiment of [37] 
(11.06 nm and 20.65 nm on the sphere and plate, respec- 
tively [37]). However, they are much lower than the highest 
peaks in the experiment [36j . The respective zero rough- 

(s) (v) 

ness levels on the sphere and on the plate, and Hq , 
are found from 



(j>,«) 



K 

£ 

»=i 



H, 



= 0. 



(8) 



From ([8]) using the roughness data one obtains Hq = 

5.01 nm and Hq — 9.66 nm. Note that precise measure- 
ments of absolute separations z discussed above result in 
separations just between the zero roughness levels deter- 
mined in (|8j). 

Special tests were performed to investigate possible 
nonlinear behavior of the oscillator under the influence 
of the Casimir force. First, the resonance frequency uj r 
observed under the excitation leading to a harmonically 
varying separation with amplitude A z was compared with 
the resonance frequency with no excitation (i.e., with sep- 
aration varied just through the thermal noise). When the 
amplitude A z was less than 4nm, no deviation was ob- 
served between the measured resonance frequency and the 
thermal resonance frequency within the rs 5 mHz noise. 
This was performed at different separations (recall that 
the amplitude actually used in the experiment was A z « 
2 nm) . The value of Af at which deviations are observed 
is a function of separation. For example, at z = 199.8 nm 
Af = 4.5 nm, at z = 247.3 nm Af = 10.0 nm, and at 
z — 302.4 nm, A"' = 15.0 nm. In all cases, when observed, 
nonlinearities decrease the resonance frequency. 

Another test performed was a check for the strength 
of the signal at different harmonics of the excitation. The 
experiment was done with the excitation z(t) at the res- 
onant frequency uj r . The checks were performed with the 
excitations at frequencies 2oj r and u> r /2, but no change 
in the response at u> r was observed for A r < 4nm. For 
larger amplitudes the results were consistent with what 
was observed in the first test. These tests all verify that 

p (fiQ cm) 



50 100 150 200 250 300 350 400 

T(K) 

Fig. 2. Resistivity of the Au films (measured with an error of 
about 2%) as a function of temperature. 
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the oscillator was in fact operating in a linear regime for 
our measurements. 

It is significant that the comparison of experimental 
data with theory of the Casimir force requires knowledge 
of the optical and electronic parameters of the Au layers. 
In previous experiments all of these parameters, including 
the plasma frequency io v and relaxation parameter 7(T), 
were taken from tables |41j . For a more conclusive com- 
parison of this experiment with different theoretical ap- 
proaches, we measured the resistivity p of the Au films as 
a function of temperature in the region from T\ — 3 K to 
400 K. These measurements were performed using a four 
probe approach on Au films of the same thickness which 
were deposited at the same time as the Au deposition on 
the oscillator, and on the same substrates. The samples 
were approximately 1 mm long and 10 /im wide. The re- 
sistivity of each sample was found by taking into account 
its geometrical factor with an error of about 2% arising 
from the errors in measuring of the sample's geometry. 
The experimental data for the resistivity versus temper- 
ature are presented in Fig. 2. These data at T » Td/4 
(where To = 165 K is the Debye temperature for Au) were 
fitted to a straight line [57] 



P(T) = 



4ir 



CT 



(9) 



Here t(T) = l(T)/v F is the relaxation time, l(T) ~ T 

1/2 

is the mean free path of an electron, vp ~ w p is the 
Fermi velocity, and C = const. The fit results in C/tUp^ 2 = 
(8.14 ± 0.16) n.!?cmK~ . On the other hand, using the 
resistivity data for pure Au as a function of tempera- 
ture 58 and the previously used value of the plasma fre- 



-3/2 _ 



.00. 



qucncy uj p = 9.0 eV |4T1|59] we obtain C/ujp 
As a result we find for the Au film used in our exper- 



iment uj p = (8.9 ± 0.1) eV. Here, the absolute error of 
0.1 eV arises from the errors of the resistivity measure- 
ments. Some of the theoretical approaches to the thermal 
Casimir force require a knowledge of the relaxation pa- 
rameter. The smooth Drude extrapolation of the imagi- 
nary part of the Au dielectric permittivity, given by the 
tabulated optical data [41] . yields the relaxation parame- 
ter at room temperature 7 = 0.0357 eV (which compares 
with 7 = 0.035 eV used in previous work [5dT[3T1[5^] 1 . 



3 Comparison of experimental data with 
different theoretical approaches to the 
thermal Casimir force 

The theoretical description of both the van der Waals and 
Casimir pressures between planar plates at temperature 
T in thermal equilibrium is given by the Lifshitz formula 



[60 



p(*)=-—y) (1-5*0) / fcjdfcxffl 

77 1=0 V 1 J Jo 
x {[r T ^(6,fc±)e 29!2 -l] _1 (10) 



[r^(6,fcj> 2 ««-l] 



Here k± = |fej_| is the magnitude of the wave vector com- 
ponent in the plane of the plates, qf = k\ + (,?/c 2 , & = 
2irkBTl/Ti are the Matsubara frequencies, I = 0, 1, 2, . .. , 
Sim is Kronecker's delta symbol, and fee is the Boltzmann 
constant. The reflection coefficients for two independent 
polarizations of the electromagnetic field (the transverse 
magnetic and transverse electric) are defined as 



?"tm (£( j k±) = £iqi f' , t T e(6: 



em + h ' 



ki+qi' 



where 



£1 







k 2 , 



mo, 



(11) 



(12) 



and e{u>) is the frequency-dependent dielectric permittiv- 
ity of the plates. 

Note that (fT0|) is the expression for a plate of infinite 
thickness. Using the Lifshitz formula for layered structures 
[47], it is easy to see that for Au layer thicknesses larger 
than 150 nm (as in our case) at, e.g., z — 400 nm the error 
due to the replacement of a layer with a semispace is less 
than 0.003%. 

It is known that there is some controversy concerning 
the application of (|10p. ITT]) to real metals. These contro- 
versies arise from different approaches to the calculation 
of the zero- frequency (I = 0) term in (jTTJ|> . For real mate- 
rials (Au for instance) e(i£i) is usually found through the 
Kramcrs-Kronig relation 



= 1 + -p 

7T 



loe"(lo) 
« 2 +& 



(13) 



where e"(w) is the imaginary part of the dielectric permit- 
tivity and the integral is taken as a principal value. Optical 
data for e"(uj) are available within some restricted fre- 
quency region [41] , and it is common to smoothly extrap- 
olate available data to lower frequencies using the imagi- 
nary part of the Drude model dielectric permittivity 



lu(lo 2 + 7 2 ) 



(14) 



If such an extrapolation is performed down to lower fre- 
quencies, including zero frequency [44] . the use of the re- 
sulting e(i£,i) in the Lifshitz theory leads to a violation of 
the Nernst heat theorem for perfect crystal lattices [61] . 
and the Casimir pressures calculated using (|10[) are in con- 
tradiction with experiment [36 ,37,39]. Because of this, two 
other approaches to the determination of e(it;i) were pro- 
posed in the literature. According to the plasma model 
approach [62]I63] . the tabulated optical data are not used 
and e(it;i) is found from the free electron plasma model 



1 + — 



(15) 



According to the impedance approach [64, 65j , the reflec- 
tion coefficients (|lip are expressed in terms of the Leon- 
tovich surface impedance Z(uj) instead of the dielectric 
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permittivity. The contributions of all Matsubara frequen- 
cies with I > 1 are obtained from tabulated optical data 
extrapolated by the Drude model using a relation Z(i^i) — 
1 / \Je{i^i). This leads to approximately the same calcula- 
tion results as the use of the dielectric permittivity. As 
to the contribution of zero Matsubara frequency it is ob- 
tained using the impedance of infrared optics, and is dif- 
ferent from that obtained using the Drude model (a dis- 
cussion of different approaches can be found in [66jl67]). 
Note that the impedance approach was used for the first 
comparison of the measurement data of this experiment 
with theory [39] . 

Although the plasma model and impedance approaches 
are in agreement with thermodynamics, neither can be 
considered as universally valid. The plasma model ap- 
proach completely neglects dissipation. Because of this, 
it is in agreement with measured data only for experi- 
ments 36 , 37 , 39J performed at separations larger than the 
plasma wavelength X p . As for the impedance approach, 
it is not applicable to short-separation experiments [30j 
because when z < X p the Leontovich impedance bound- 
ary conditions become invalid due to the violation of the 
inequality |Z(o;)| -C 1. 

Recently [40] a new approach to the thermal Casimir 
force between real metals was proposed which is equally 
applicable at both small and large separations. This ap- 
proach is based on the use of the generalized plasma-like 
dielectric permittivity 

eH = l-4 + E^ » (16) 

uj 2 ^— ' to 2 -uj 2 - iqauj 
j=i j *j 

which takes into account the interband transitions of core 
electrons. Here ujj ^ are the resonant frequencies of the 
core electrons, gj are the respective relaxation frequen- 
cies, fj are the oscillator strengths, and K is the number 
of oscillators. Note that the term —uj 2 /uj 2 on the right- 
hand side of ([T6j) describes free electrons and leads to a 
purely imaginary current. This contribution to e(u)) is en- 
tirely real and does not include dissipation. Importantly, 
the oscillator term on the right-hand side of (fl6|) does not 
include the oscillator with zero resonant frequency loq = 0, 
which is equivalent to the Drude dielectric function, i.e., it 
does not describe conduction electrons but only core elec- 
trons. This term incorporates dissipation due to interband 
transitions. 

In [JD] the Lifshitz theory together with the dielec- 
tric permittivity (fTl)]) was used to calculate the thermal 
Casimir force in a short-separation experiment [30], and 
the experimental results were found to be in good agree- 
ment with theory. For this purpose the oscillator param- 
eters of Au in (|T6]) were taken from [421143] where they 
were found using the 3-oscillator model fitted to old DESY 
data. Below we compare the 3-oscillator fit of [32TPI3"] with 
the more complete data set of [41] and perform a more ex- 
act 6-oscillator fit. The resulting oscillator parameters are 
used to calculate the Casimir pressure in the most precise 
experiment described in the previous section. 



The Kramers-Kronig relation ([T3"]) was derived [55] for 
dielectric permittivities which were regular or which had 
a first order pole at zero frequency. For the dielectric per- 
mittivity P^]) which has a second order pole at ui = 0, 
the Kramers-Kronig relation expressing e(iuj) in terms of 
e"(uj) is the following [40] : 

In the tables of [U] the most complete data are col- 
lected for real, ni(uj), and imaginary, 712(0;), parts of the 
complex refraction index of Au in the frequency region 
from0.125eVto 9919eV (leV = 1.519xl0 15 rad/s). From 
these data, the imaginary part of the Au dielectric permit- 
tivity is expressed as 2ni(uj)n2(uj) . To obtain the contri- 
bution of core electrons to the dielectric permittivity, we 
consider the difference 

~ 2 ~ 

e'M = 2rciMn 2 M - ^ ( J^- 2) , (18) 

where in accordance with (|14[) the subtracted term ap- 
proximately describes the contribution of free conduction 
electrons to optical data. In Fig.[3]the quantity e" Au is plot- 
ted as a function of w within the frequency region from 
2.0 eV to 25 eV (solid line). For uj < 2 eV the dielectric per- 
mittivity is determined by free conduction electrons, and 
for uj > 2.5 eV there is already practically no contribution 
from conduction electrons, and e Au (uj) xt 2ni(ui)n2(uj). 
The upper limit of the region under consideration is de- 
termined by the frequencies contributing to the Casimir 
pressure (|10|) . Even at the shortest separation considered, 
z = 160 nm, the characteristic frequency Q c = c/(2a) « 
0.62 eV. Bearing in mind that even for precise computa- 
tions of the pressure it is sufficient to take into account 
the contribution from Matsubara frequencies up to 15J? C , 
setting the upper limit of our region equal to 25 eV is more 
than adequate. 




5 10 15 20 25 



w(eV) 

Fig. 3. Tabulated optical data for the imaginary part of Au di- 
electric permittivity [41] (with the contribution of conduction 
electrons subtracted) are shown by the solid line. The oscilla- 
tor fits are shown as long-dashed line [421143] (DESY data, 3 
oscillators) and as short-dashed line (6 oscillators). 
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The solid line in Fig. [3] was fitted to the imaginary part 
of dielectric permittivity IT6|) 

e» = y . a 2 ? (19) 

with X = 6 oscillators. The resulting set of oscillator pa- 
rameters fj, LUj and gj is presented in Table 2. In Fig.[5]thc 
imaginary part of permittivity calculated using the ana- 
lytic expression (fT9]l is shown by the short-dashed line. In 
the same figure the 3-oscillator fit from [42ll43j is shown 
as the long-dashed line. As is seen in Fig. [3] the short- 
dashed line based on the 6-oscillator fit better reproduces 
the actual data than does the long-dashed line using the 
3-oscillator fit. 

Table 2. The oscillator parameters for Au in equations (|16[1 
and (|19[) found here from the 6-oscillator fit to the tabulated 
optical data in |41] , 



3 


Ui (eV) 


95 (eV) 


fi (eV 2 ) 


1 


3.05 


0.75 


7.091 


2 


4.15 


1.85 


41.46 


3 


5.4 


1.0 


2.700 


4 


8.5 


7.0 


154.7 


5 


13.5 


6.0 


44.55 


6 


21.5 


9.0 


309.6 



The Casimir pressure Pl(z) at all separations of inter- 
est was computed using the Lifshitz theory in equations 
([T0" ]) -([l2" ]) and (16]) with uj p = 8.9 eV, as determined for 
our films in Sec. 2, and the oscillator parameters from 
Table 2. For comparison with the experimental data, the 
values of Pl(z) were geometrically averaged over all possi- 
ble separations between the rough surfaces weighted with 
the fractions of the total area occupied by each separa- 
tion, as discussed in Sec. 2. This results in the theoretical 
Casimir pressures taking surface roughness into account 
via the equation 

p th &) = E E <4S W 
k=i j=i 

x P L [zi + + H { P) - 4 S) - hf) . (20) 

The pressures P th (zi) were computed at each experimen- 
tal point Zi. Note that (|20|) takes into account the com- 
bined (nonmultiplicative) effect of nonzero temperature 
and finite conductivity on the one hand [this is incorpo- 
rated in Pl(z) computed using the Lifshitz formula], and 
of surface roughness on the other. The contributions of 
diffraction-type and correlation effects in the roughness 
correction [691170] , which are not taken into account in the 
geometrical averaging (|2H]) . were shown to be negligible 
[37] . In this experiment the contribution of the roughness 



correction to the Casimir pressure computed using (|20|) 
is very small. For example, at z = 162 nm the roughness 
correction contributes only 0.52% of the total pressure. 
At separations z = 170, 200 and 350 nm roughness con- 
tributes only 0.48, 0.35 and 0.13% of the Casimir pressure, 
respectively. The magnitudes of the computed theoretical 
Casimir pressures at some experimental separations are 
listed in column (b) of Table 1. 

We now discuss the accuracy of our computations. One 
of the sources of the theoretical errors is the sample-to- 
sample variation of the optical data for the complex in- 
dex of refraction. As was shown in [37] (see also [3pJ), in 
our experiments the variation of the optical data leads to 
an uncertainty in the magnitude of the Casimir pressure 
which is substantially smaller than 0.5%. To be conserva- 
tive, we admit an uncertainty as large as 0.5% in the com- 
putations due to the use of tabulated optical data over the 
entire measurement range. There are claims in the litera- 
ture [71] that the theoretical computations of the Casimir 
pressure between gold surfaces are burdened by up to 5% 
errors due to the use of different Drude parameters mea- 
sured and calculated for different samples. This is, how- 
ever, irrelevant to our experiment. The hypothesis that the 
magnitude of lu p is much smaller than the value we have 
used above (i.e., uj p = 6.85 eV or 7.50 eV, as suggested in 
[71 j ) is rejected at high confidence by our experiment, and 
by all previously performed measurements of the Casimir 
force between Au surfaces. 

The other possible source of theoretical errors is con- 
nected with the fact that we compute the Casimir pressure 
at experimental separations which are determined with an 
error Az = 0.6 nm [72]. Noting that the dominant the- 
oretical dependence of the Casimir pressure is z~ 4 , one 
finds that the relative error in the pressure is equal to 
AAz/z. This varies from 1.5% at z = 160 nm to 0.32% 
at z = 750 nm. The other theoretical errors, e.g., aris- 
ing from neglect of patch potentials or spatial nonlocality, 
were analyzed in detail in [37] and found to be negligible. 
By combining the above two theoretical errors discussed 
here and in the previous paragraph at a 95% confidence 
level using the statistical procedure applicable to system- 
atic errors described by a uniform distribution 37,48], we 
obtain the total theoretical error A tot P th (z) as a func- 
tion of separation. This error assumes a maximum value 
of 18.7 mPa at z — 162 nm, which is almost 9 times larger 
than the total experimental error. Note that in a simi- 
lar analysis in [37j . one additional theoretical error due 
to the use of PFA was included. In [37] it was first com- 
bined with the theoretical error due to sample-to-sample 
variation of the optical data, with the result that the dis- 
tribution law of the combined quantity was not uniform. 
In this work, however, the error due to the accuracy of 
PFA is included with the experimental errors. Because of 
this, the total theoretical error is determined by only two 
contributions. With the increase of separation to z = 300, 
400 and 746 nm, the total theoretical error decreases to 
1.15, 0.34 and 0.024mPa, respectively. The relative theo- 
retical error AP th (z)/\P th (z)\ assumes a maximum value 
of 1.7% at z = 162 nm. When the separation increases 
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to z — 300, 400 and 746 nm, the relative theoretical er- 
ror decreases to 1.0, 0.86 and 0.65%, respectively. This is 
mainly explained by the decreased role of uncertainty in 
determining the separations. 

We can now compare experiment and theory by con- 
sidering the differences P th (zi) — P{zi) at each experimen- 
tal separation Z{. The confidence interval for the quantity 
P th (zi) — P(zi) determined at 95% confidence probability 
is given by [— So. 95(24), So. 95(24)] where the half-width of 
this interval can be found using the composition rule [37|, 
EE] 

So.w(*) = min {A tot P th ( Zl ) + A tot P e ^( Zi ), 

fcoW [A tot P ih (zi)] 2 + [A**P°*»(zi)] 2 \ . (21) 

(2) 

Here for two composed quantities fc 95 = 11- The values 
of the half-width of the confidence interval are listed in 
the last column of Table 1. 

In Fig. 4a the differences P th (zi) — P(zi) at all exper- 
imental points are shown as dots. In the same figure the 
confidence interval [—So. 95(24), So. 95(24)] at each z is sit- 
uated between the solid lines. As seen in the figure, all 
dots (and not only 95% of them as required by the rules 
of mathematical statistics) are well inside the confidence 
interval at all separations considered. This means that the 
experimental data are consistent with theory based on the 
generalized plasma- like dielectric permittivity (|16p . and 
that in our conservative error analysis the errors (espe- 
cially at short separations) are overestimated. For com- 
parison purposes in Fig. 4b we plot as dots the differences 
P th (zi) — P(zi) where the experimental data are the same 
as in Fig. 4a, but with P th (zi) computed as in [35] us- 
ing the Leontovich surface impedance approach, with the 
Drude parameters lo p = 8.9 eV, 7 = 0.0357 eV. In Table 1, 
column (c) we present the magnitudes of the Casimir pres- 
sures P th (zi) computed using the surface impedance ap- 
proach at different separations. As is seen in Fig. 4b, the 
impedance theoretical approach is also consistent with the 
data. However, while in Fig. 4a there are practically no de- 
viations between experiment and theory at z > 350 nm, 
in Fig. 4b the deviations are noticeable up to z = 450 nm. 
By comparing columns (b) and (c) in Table 1, we can 
conclude that the differences between the two theoretical 
approaches do not exceed the magnitude of the theoretical 
error. The comparison of columns (6) and (c) with column 
(a) shows that at all separations the approach using the 
generalized plasma-like model is in somewhat better agree- 
ment with data than the surface impedance approach. As 
is seen in Fig. 4a, b, the largest deviations between both 
theoretical approaches and experimental data are at short 
separations from 162 to 200 nm. Although these deviations 
are not statistically meaningful, because they are well in- 
side the confidence interval, in Sec. 5 we will discuss possi- 
ble reasons leading to the deviations between experiment 
and theory at shortest separations. 

A completely different situation occurs when we com- 
pare the experimental data with the alternative approach 
to the thermal Casimir force |44j using the Drude model 



ptn _ p ( mPa ) 
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Fig. 4. Differences between theoretical Casimir pressures com- 
puted using the generalized plasma model approach (a) and the 
Leontovich surface impedance approach (b) and mean experi- 
mental Casimir pressures (dots) versus separation. Solid lines 
indicate the limits of the 95% confidence intervals. 



to compute the contribution of the zero- frequency term in 
the Lifshitz formula. To perform the comparison, we cal- 
culate the theoretical Casimir pressures in the framework 
of [34] with the refined values of the Drude parameters 
Up = 8.9 eV, 7 = 0.0357 eV (all details of this approach 
and of computations can be found in [37]). The magni- 
tudes of the resulting Casimir pressures P^ at a few differ- 
ent separations are listed in Table 1, column (d). In Fig. 5 
we plot the differences Pj}(zi) — P(zi) at all experimental 
separations. The confidence interval [—So. 95(24), So. 95(24)] 
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at each zi is the same for all theoretical approaches. Once 
again, the limits of the confidence interval are denoted by 
the solid lines (in Fig. 5 only one solid line is shown be- 
cause practically all dots are above it). As is seen in Fig. 5, 
the Drude model theoretical approach is experimentally 
excluded at a 95% confidence level within the whole mea- 
surement range from 162 to 746 nm. This conclusion is 
confirmed by the calculation data in Table 1 . Subtracting 
the magnitudes of the theoretical Casimir pressures, \P%\, 
in column (d) from the experimental results, |P|, in col- 
umn (a), we obtain at all separations larger results than 
the half- width of the confidence interval, £0.95 (z), given 
in column (e). 

The wide gaps between the solid line and dots in Fig. 5 
suggest that the Drude model approach is actually ex- 
cluded experimentally at an even higher confidence than 
95%. To make this argument quantitative, we calculate 
the half-width of a confidence interval at 99.9% confidence 
from 

£b.999(2) _ £(l+0.999)/2 

(32) 



■=0.95(2) t(i+Q.95)/2(32) 



1.85, 



(22) 



where t p (f) is the Student coefficient used in Sec. 2. The 
limits of the 99.9% confidence intervals obtained in (|22l) 
are shown in Fig. 5 by the dashed line. As is seen in Fig. 5, 
the differences P^ — P are found outside of the 99.9% con- 
fidence interval at separations from 210 to 620 nm. This 
conclusively demonstrates that our experiment is irrec- 
oncilable with the Drude model approach to the thermal 
Casimir force. At the same time, the approaches based on 
the generalized plasma-like dielectric permittivity, and on 
the Lcontovich surface impedance, are consistent with ex- 
periment. Note that in our experiment the Drude model 
approach is excluded at separations below 1 /xm. In the 
proposed experiments [73"II74"U75"] it is planned to test the 
predictions of different theoretical approaches to the ther- 



pth 
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Fig. 5. Differences between theoretical Casimir pressures com- 
puted using the Drude model approach and mean experimental 
Casimir pressures (dots) versus separation. The solid line in- 
dicates the limits of the 95% confidence intervals while the 
dashed line indicates the limits of the 99.9% confidence inter- 
vals. 



mal Casimir force at separations of about several microm- 
eters. 



4 Constraints on Yukawa-type hypothetical 
interactions and light elementary particles 

As was mentioned in the Introduction, at separations be- 
tween macroscopic bodies of about 1 /im and less, the Ca- 
simir force is the dominant background force. From the 
level of agreement between the experimental data for the 
Casimir pressure and Lifshitz theory (with a generalized 
plasma- like permittivity in Sec. 3), one can constrain any 
additional force which may coexist with the Casimir force. 
As noted in the Introduction, many extensions of the stan- 
dard model predict a Yukawa correction to the Newtonian 
potential energy between two point masses mi and ui2 at 
a separation r, given by [TUllT^lfTS] 



V(r) = -=—=■ (1 + ae r/A J 



(23) 



Here G is the Newtonian gravitational constant, a is a 
dimcnsionless constant characterizing the strength of the 
Yukawa interaction, and A is its range. 

The total force acting between two parallel plates due 
to the potential (|2"3"|) can be obtained by integration of (|2"3")l 
over the volumes of the plates, and subsequent negative 
differentiation with respect to z. In experiments measur- 
ing the Casimir force the contribution of the gravitational 
force is very small and can be neglected 35 , 36J . Thus, in 
what follows we consider only the contribution from the 
Yukawa term in (|2"3")l . 

To find the Yukawa pressure for our setup we should 
take into account the detailed structure of our test bod- 
ies. (As was shown in Sec. 3, for the calculation of the 
Casimir pressure it is possible to replace the Au coating 
films with Au semispaces and we need not consider the un- 
derlying substrate.) In the present experiment a sapphire 
sphere of density p s = 4.1g/cm was first coated with a 
layer of Cr of density p c — 7.14g/cm and thickness A c = 
10 nm, and then with an external layer of gold of thickness 
Ag } = 180 nm and density p g = 19.28 g/cm 3 . The Si plate 
of thickness L — 3.5 /im, and density psi = 2.33 g/cm 3 
was also first coated with a layer of Cr of A c = 10 nm 
thickness, and then with a layer of gold of A g = 210 nm 
thickness. Under the conditions z, A <C R, satisfied in this 
experiment, the equivalent Yukawa pressure between the 
two parallel plates with the same layer structure as the 
above sphere and a plate is given by [32 . 47 



P h yp( z ) = -2TrGa\ 2 e~ z / x (24) 
x p. ( Pg Pc )e-^" - ( Pc Pa )e-^^y\ 
x 'p g (p g - Pc ) e -4 P, /A _ {pc _ ps . )e -(^ + A c yx] 



We have verified that the surface roughness, as reported in 
Sec. 2, cannot significantly affect the magnitude of a hy- 
pothetical pressure with an interaction range longer than 
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10 rim, and hence can be neglected. Because of this, there 
is no need to perform geometrical averaging as in (|20| 
when calculating the Yukawa interaction. 

According to Sec. 3, theories of the thermal Casimir 
force using the generalized plasma-like permittivity or the 
Lcontovich surface impedance are consistent with exper- 
imental data. As was noted in Sec. 3, in our conserva- 
tive analysis the errors [and consequently the width of 
the confidence interval 2S(z)] are overestimated. The rea- 
son for this is that we have included the error due to the 
uncertainty of experimental separations Az in the anal- 
ysis of the theoretical errors. As a result, the theoretical 
pressures acquired an extra error of kAAz/ z which led 
to enormous widening of the confidence interval at short 
separations (see Figs. 4 and 5). This approach was use- 
ful in selecting among different theories of the thermal 
Casimir force, and permitted us to exclude the one based 
on the use of the Drude model at practically 100% con- 
fidence. However, as is clearly seen in Fig. 4a, b, the ac- 
tual width of the confidence interval is much less than 
that between the solid lines (recall that the actual confi- 
dence interval determined at 95% confidence should con- 
tain about 95% of the data dots but not all of them). It is 
easily seen that if the theoretical error 4Az / z due to un- 
certainties in experimental separations is disregarded, the 

resulting more narrow confidence interval —S{z),E(z) 

still contains all dots representing P th (z) — P(z) within 
the separation region from 180 to 746 nm. At a separa- 
tion z — 180 nm, the half-width E = 4.80 mPa. At typical 
separations z = 200, 250, 300, 350, 400 and 450 nm S is 
equal to 3.30, 1.52, 0.84, 0.57, 0.45 and 0.40 mPa, respec- 
tively. Thus, for 180 nm < z < 746 nm the magnitude of 
the hypothetical pressure should satisfy the inequality 



\P hyp (z)\ < E(z). 



(25) 



Bearing in mind that the half- width of the confidence in- 
terval E(z) was defined at a 95% confidence, the same 
confidence also applies to the constraints following from 
the inequality (|25|) . 

We have performed a numerical analysis of equations 
(J2U) and (p2"5)) at different separations and determined the 
resulting region of (A, a)-plane where the inequality (|25j) 
is satisfied, so that the existence of Yukawa interaction 
is consistent with the level of agreement achieved between 
data on the measurement of the Casimir force and relevant 
theory. The strongest constraints within the interaction 
region 10 nm < A < 56 nm are obtained from the com- 
parison of measurements with theory at a separation z = 
180 nm. With the increase of A, the strongest constraints 
on a were obtained from the consideration of larger sep- 
arations. Thus constraints in the regions 56 nm < A < 
71 nm, 71 nm < A < 89 nm, 89 nm < A < 140 nm, 140 nm < 
A < 220 nm and 220 nm < A < 500 nm were obtained 
from the agreement between Casimir pressure measure- 
ments and theory at separations z = 200, 250, 300, 350 
and 400 nm, respectively. 

The resulting constraints on a are plotted in Fig. 6 for 
different values of A (line 1). The region in the (A, a)-plane 
above the line 1 is excluded by the results of the Casimir 



pressure measurements compared with theory, and below 
line 1 is allowed. For comparison, constraints from earlier 
experiments are also shown in Fig. 6 in a similar manner. 

Special attention should be paid to line 2 represent- 
ing constraints following from the short-separation exper- 
iment [30] on the measurement of the Casimir force be- 
tween an Au-coated sphere and a plate using an atomic 
force microscope. The constraints on a Yukawa hypothet- 
ical interaction following from that experiment were ob- 
tained in [35], and later used in [36 p 37 | l39] for comparison 
with constraints following from other experiments. How- 
ever, in [35j the level of agreement between experiment 
and theory at zero temperature was described in terms of 
the root-mean-square deviation which, as was recognized 
later in [37], is not an appropriate quantity in strongly 
nonlinear situations. In addition, the calculational scheme 
using the root-mean-square deviation does not permit us 
to determine the confidence level of the results. 

Here we reanalyze the experimental data of [30j and 
compare them with theory using the Lifshitz formula at 
laboratory temperatures in a sphere-plate configuration, 
supplemented by the generalized plasma-like dielectric per- 
mittivity (|16p . The results of this reanalysis are expressed 
in terms of the confidence interval [— O(z), 0{z)] deter- 
mined at 95% confidence for the differences between the- 
oretical and mean experimental Casimir forces, F (z) — 
F(z). This interval includes, in particular, the theoretical 
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Fig. 6. Constraints on the strength of the Yukawa interaction 
versus interaction range. Line 1 is obtained in this paper, line 
2 was obtained in [35] using the Casimir force measurement of 
[30] and adapted in this paper to the accepted 95% confidence 
level. Lines 3 and 4 were obtained in [78] and [37], respectively. 
Line 5 was obtained in the first reference of 32 using the 
Casimir force measurement [27]. The region of the (A, a) plane 
above each line is excluded and below the line is allowed. 
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errors 3Az / z arising due to uncertainties of experimental 
separations in the sphere-plate configuration. This inter- 
val cannot be narrowed as we did above in the case of the 
present experiment, because the measurement in |30j is 
inherently noisier. For example, at z = 61.08 nm the half- 
width of the confidence interval is — 31.6 pN, and with 
an increase of separation up to 100.15 and 200.46 nm it 
decreases to 9.17 and 7.20 pN, respectively. The resulting 
constraints at a 95% confidence level are determined from 

|F h ^(z)| <G(z), (26) 

where F h yP(z) is the Yukawa hypothetical force acting 
between an Au coated sphere and a plate [3"5"lf3"rJ] . These 
constraints are represented by line 2 in Fig. 6. Note that 
the constraints given by line 2 are up to order of mag- 
nitude weaker than those in [35], but they benefit from 
high confidence, and can be compared with future work 
on the subject by using the same rigorous approach to the 
comparison experiment with theory as proposed in [37, 76, 
[77]. 

The other lines in Fig. 6 are obtained from the Casimir- 
less experiment [78] (line 3), previous measurements of the 
Casimir pressure using the micromachined oscillator [37j 
(line 4), and in [32] (the first paper) from a torsion pendu- 
lum experiment [2?] (line 5). As is seen in Fig. 6, the result- 
ing constraints represented by line 1 are strongest within 
the interaction range 20 nm < A < 86 nm with the largest 
improvement by a factor 4.4 at 26 nm. Note that further 
strengthening of the resulting constraints on a within a 
submicrometer interaction range could provide important 
information concerning predicted particles such as scalar 
axions, graviphotons, hyperphotons, dilatons, and moduli 
among others. For such particles the interaction constant 
a could be much larger than unity. The same holds for the- 
ories based on extra-dimensional physics with low-energy 
compactification scale where, for instance, for models with 
three extra dimensions the predicted characteristic size of 
extra dimensions is of about 5 nm [T2HT3] . 

To conclude this section we briefly discuss possible 
reasons for the observed deviations between experiment 
and theory at the shortest separations shown in Fig. 4a,b. 
These deviations are well inside the 95% confidence inter- 
val determined for P th — P and thus they are not statis- 
tically meaningful. Nevertheless if we bear in mind that 
the deviations under consideration are several times larger 
than the total experimental error, there may be some 
underlying physics leading to the small discrepancies be- 
tween experiment and theory. The most natural assump- 
tion is that there is some undiscovered nonlinearity of the 
oscillator which results in an additional systematic error at 
short separations. However, as discussed in Sec. 2, special 
tests of the oscillator linearity have been performed which 
did not reveal a nonlinear behavior for the amplitudes of 
sphere oscillations employed in this experiment. Another 
possible effect may be connected with some fine properties 
of interacting surfaces determined, e.g., by correlation ef- 
fects in surface roughness or by patch potentials. However, 
as was analyzed in detail in [37] , these effects are negligibly 
small. Thus, even assuming enormously large patches due 



to monocrystals with grain sizes ranging from ~ 300 nm 
(i.e., larger than the film thickness) and to 25 nm, the cor- 
rection to the pressure due to patches at z = 160 nm is 
only 0.42 mPa (to be compared with the deviation between 
experiment and theory of almost lOmPa, as in Fig. 4a). 

We next consider the possibility that the deviation 
may be caused by the Yukawa interaction (|23p with some 
appropriate values of a and A. A simple calculation shows 
that the deviations between experiment and theory at 
short separations would practically disappear if we allowed 
a Yukawa interaction with a = 5.0 x 10 21 and A = 10.4 nm. 
This interaction would correspond to a point in Fig. 6 
with logA = —7.98 situated slightly below line 1. How- 
ever, such a point would lie above the point a = 1.0 x 10 20 
with the same A on line 2, which implies that the assumed 
Yukawa interaction is excluded by the AFM experiment 
[30j . Bearing in mind that in the above we have rean- 
alyzed the results of [30] at a 95% confidence level us- 
ing modern methods of comparison between experiment 
and theory, this is strong evidence against the existence 
of a single Yukawa interaction with a = 5.0 x 10 21 and 
A = 10.4 nm. However, this analysis does not necessarily 
exclude the possible existence of more than one Yukawa, or 
other interactions having a different spatial dependence. 
A more decisive conclusion about the presence of hypo- 
thetical interactions can be obtained through a repetition 
of the experiment described in Sec. 2 using a Si plate but 
with no covering metallic layer. Noting that the density 
of Au is in 8.3 times larger than the density of Si, and 
that the Casimir pressure between Au and Si is approx- 
imately 1.5 times smaller than between Au and Au, the 
deviation caused by the Yukawa interaction should prac- 
tically disappear if the Au coated plate is replaced with a 
Si plate. 

5 Conclusions and discussion 

In this paper we have presented additional details on the 
recent experimental determination of the Casimir pres- 
sure between two parallel plates using a micromachined 
oscillator. This experiment incorporates several improve- 
ments over all previous measurements. In particular, the 
measurements over a wide separation range were repeated 
many times at practically the same points for each repeti- 
tion. This permits us to substantially reduce the random 
error, and to make it much smaller than the systematic er- 
ror for the first time in Casimir force measurements. Also 
the plasma frequency of the Au films was determined using 
the measured temperature dependence of their resistivity. 

The resulting experimental data were compared with a 
new theoretical approach to the thermal Casimir force us- 
ing the Lifshitz theory incorporating a generalized plasma- 
like dielectric permittivity which takes into account the 
intcrband transitions. For this purpose a new oscillator 
fit of the tabulated optical data for the imaginary part of 
the dielectric permittivity of Au was performed which is 
more exact than a previously used fit based on DESY data. 
The new theoretical approach was also compared with the 
previously known approach using the Leontovich surface 
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impedance, and with the alternative approach using the 
Drude model. The Drude model approach was excluded 
experimentally at a 99.9% confidence level over a wide 
separation range. 

One of the main aims of this paper is the applica- 
tion of the Casimir effect to obtain stronger constraints 
on hypothetical long-range interactions and light elemen- 
tary particles. We have reanalyzed the previously known 
constraints from the measurement of the Casimir force 
between an Au-coated sphere and a plate using modern 
methods of comparison of experiment and theory at high 
confidence. We have also used the resulting level of agree- 
ment between the measurements of the Casimir pressure 
and the new theory to strengthen constraints on the hy- 
pothetical Yukawa-type interaction. The new constraints 
obtained above are the strongest within the interaction 
range from 20 to 86 nm, with the largest improvement by 
a factor 4.4. These results are relevant for the verification 
of different theoretical predictions made on the basis of 
unified field theories beyond the standard model, and of 
extra-dimensional physics. We have also discussed some 
possible reasons for small deviations between experiment 
and theory at the shortest separations considered. It was 
shown that although these systematic deviations are not 
statistically significant, the fact that we have no expla- 
nation for them at present suggests that further experi- 
mental and theoretical work is required to elucidate their 
nature. 
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